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ABSTRACT 

Cold clouds embedded in warm media are very common objects in astrophysics. Their 
disruption timescale depends strongly on the dynamical configuration. We discuss the 
evolution of an initially homogeneous cold cloud embedded in warm turbulent gas. 
Within a couple of dynamical timcscalcs, the filling factor of the cold gas within the 
original cloud radius drops below 50%. Turbulent diffusivities estimated from the time 
evolution of radial filling factor profiles are not constant with time. Cold and warm 
gas are bodily transported by turbulence and mixed. This is only mildly indicated 
by column density maps. The radiation field within the cloud, however, increases by 
several orders of magnitudes due to the mixing, with possible consequences for cloud 
chemistry and evolution within a few dynamical timescales. 

Key words: hydrodynamics — radiative transfer — turbulence — ISM: clouds — 
ISM: kinematics and dynamics. 



1 MOTIVATION 

Molecular clouds (MCs) in the Galaxy exhibit a wealth 
of structures in (column) densities, velocities and mag- 
netic fields. The observed structural properties strongly 
suggest that MCs are highly dynamical objects within a 
turbulent interstellar medium (ISM). Especially, the ob- 
served non-thermal linewidths iFalgarone fc Philipsl Il990t 
IWilliams et~aTl l2000tl indicate the turbulent nature of the 
clouds. The importance of molecular cloud turbulence 
for the process of star for mation has been the sub ject 
of many investigations ( see lElmegreen fc Scald 120041 and 
iMac Low fc K lcssen 20Q4| for overviews). 

Turbulent mixing is an ever-recurring theme in the ISM. 
Classical turbulent mixing accelerates diffusive transport by 
a fa ctor on the order of the Reynolds numbe r of the flow 
fe.g ISchatzmarJll977l).lde Avillez fc Mac Low! \200l 12003ft 
and lKkssen^^AiT(l2003l) discussed the applicability of tur- 
bulent transport and mixing to a turbulent ISM, stating 
that transport rates can vary strongly with time, some- 
times exhibiting super-diffusive behaviour due to bulk mo- 
tions of the gas. Turbulent mixing has also been held re- 
sponsible as a source of the highly ionized gas observed 
around High V elocity Clouds in the Galactic Halo (e.g. 
iFox et al.ll200^ . although it is unclear whether turbulent 
mixing or evaporation by heat conduction is the dominat- 
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ing process (ICowie fe McKedll977t iMcKee fc Cowid Il977l : 
iBalbus fe McKedll982HLazarianll200g) ~ 

A spherical cold cloud travelling through a warm uni- 
form medium dissolves approximately within a timescale 
Td = M c i/(pTrRciv), where M c i is the cloud's initi al mass, 
p its density, R c i its radius, and v its velocity iNulsenl 
Il982fl . The main agent is the Kelvin-Helmholtz- instability, 
whose efficiency in dissolving the cloud depends on the 
cloud's physical prop erties, e.g whether it is self-gravitating 
iMurrav et aljlT993l) - in which case there exists a critical 
mass above which the cloud r emains stable -, or whether it 
suffers radiative energy losses llVietri et al]ll997ll - in which 
case the instability can be damped, stabilizing the cloud. 

Cloud dispersal by hydrodynamical instabilities might 
be of considerable importance: observations and theoreti- 
cal considerations suggest that molecular clouds in the so- 
lar neighborhood have lifetimes of approximately 2 to 3 
Myrs, which would necessitate close to i nstantaneous star 
formation once the clouds have formed (l ElmeeTeen| [l993|, 
| 2000|: lHartmann et alJl200lt IPringle et al.ll200lt iHartmannl 
2003). While stellar feedback in the form of supernovae 
could disrupt a cloud, low-mass stars in solar neighbor- 
hood clouds might not be efficient enough to achieve 
such a rapid dispersal (e.g. IMac Low! I2004T) . An alter- 
native route has been offered by interpreting molecular 
clouds as transient objects generated by l arge-scale c ollid- 
ing fl ows (iBalk^teros^^ |Hjrtmann_e^_alJ 
| 200lt iHeitsch et al l 120051: IVazauez-Semadeni et al.l l2006t 
lH^tech^t^ni200c3) '! In this picture, the flows in which the 
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clouds form eventually might lead to their dispersal within 
a few Myrs. 

We investigate the role of turbulence for the overall evo- 
lution of a cold cloud - corresponding to the Cold Neutral 
Medium (CNM) - embedded in warm gas whose parame- 
ters are representative of the Warm Ionized Medium (WIM). 
Specifically, we are interested in the timescales on which a 
cold solid cloud will disperse in a turbulent environment, 
and in the evolution of the radiation field within the cloud 
O- A detailed, time-independent study of the radiation 
field i n a f ilamentary cloud ca n be found in iBethell et alJ 
J2004D and lBethell et alJ J2006I) . 

We find (|3J that the warm and cold components are 
efficiently mixed by bodily transport. However, the central 
optical depth stays comfortably above 1 even for late times. 
At first, this might seem contradictory, however, the local 
radiation field within the cloud increases by several orders 
of magnitude within a couple of dynamical times. Turbulent 
diffusivities derived from the expansion of the cold cloud 
generally are not constant with time. Possible consequences 
for cloud lifetimes and evolution are discussed in |1] 



2 THE PROBLEM AND ITS SETUP 

We start with a uniform sphere of cold dense gas in a warm 
diffuse ambient medium. The system is initially completely 
in thermal and almost in turbulent pressure 1 equilibrium. 
One model set (sequence A) is purely adiabatic, i.e. the two- 
phase structure of the system will get erased with time as 
gas with temperature intermediate to the two initial tem- 
peratures is created. The second set (sequence C) keeps the 
identities of the warm and cold phases distinct by a combi- 
nation of heating and cooling processes iWolfire et al.lll995T) 
typical for the warm and cold atomic ISM. Thus, we study 
turbulent mixing under conditions with and without radia- 
tive losses. Initially, both gas components are in thermal 
equilibrium at temperatures that for the radiative case cor- 
respond to the two stable temperature regimes, i.e. there is 
no gas in the thermally unstable regime initially. The lower 
initial temperature is set to T — 31K, and the density con- 
trast between warm and cold gas is 300, with no = 0.5cm -3 
in the warm phase at T = 9.2 x 10 3 K. The cubic box is 
periodic in all directions, with a side length, L, of 44pc, and 
the cold gas sphere starts out with a radius of 4.4pc. 

Instead of studying the cloud evolution within a shear 
flow, we resort to a more direct way of treating turbulence, 
namely by initially imposing a velocity fi eld drawn from 
a ran dom Gaussian distribution (see e.g. iMac Low et alJ 
1998). Power is alotted in Fourier space only to the largest 
scales, and with random phases. This is meant to mimic 
the effect of turbulence generated by an unspecified source 
on larger scales. The initial Mach number in the warm gas 

1 As discussed by e.g. Ballcstcros-Parcdcs (2006), the concept of 
turbulent pressure is only applicable if there exists a scale sepa- 
ration between the (large) object scale and the (small) turbulent 
scale. In a medium without scale separation - e.g. with a self- 
similar turbulent power spectrum - the term "turbulent pressure" 
should be read as "turbulent kinetic energy density" . It is in this 
sense we will use the term "turbulent pressure" throughout the 
paper. 



A4o = 2 or 3 - corresponding to vo — 22.4km s~ and 
vo = 33.6km s _1 - is higher than Mach numbers commonly 
ascribed to the warm ISM, however, since the turbulence 
is not driven, the system acquires reasonable values of M 
once turbulence is fully developed. The justification for this 
will be discussed in H3.4I One might argue that the sce- 
nario of a uniform, spherical cold cloud evolving in a decay- 
ing turbulent velocity field is only of limited physical rel- 
evance, since the Galactic ISM turbule nce is to a large ex- 
tent thought to be driven by supernovae iCoxJL-Sjnith|l974]; 
[ McKee fe OstrikeJl97'AlRosen fe BregmarJl995l : lde Avillezl 
120001 Ide Avillez fc Berrvll200ll) . However, the goal of this 
study is not to model the evolution of a cold (possibly molec- 
ular) cloud, but to demonstrate that even under unfavor- 
able conditions (decaying turbulence, no stellar feedback, no 
gravitational fragmentation, starting with a uniform spher- 
ical cloud), cold clouds fragment and disperse sufficiently 
within a dynamical timescale so that the internal struc- 
ture of the cloud is thoroughly altered. Driven turbulence 
and/or stellar feedback obviously would lead to faster dis- 
persal, while gravitational fragmentation would result in a 
smaller gas filling factor, thus opening holes and channels 
for radiation and energy (pos sibly in form of waves, see e.g. 
iHennebelle fc Inutsukal BOOG) to enter the cloud. For this 
reason also a non-uniform cloud would be more prone to 
disperse in a turbulent environment. 

The adiabatic model sequqence we will denote by A, the 
radiatively cooled model by C. The initial Mach number 
of the model is denoted by "2" or "3", for Mach 2 or 3. 
We ran models at linear resolutions of 128 3 and 256 3 cells, 
indicated by the letters a and b respectively. The choice of 
the - rather small - initial cloud radius could raise concerns 
about how well even the 256 3 -models are resolved. Thus, we 
ran an additional model, A21, with the same parameters as 
the adiabatic model at Mach 2, A2b, except for the initial 
cloud radius, here set to 8.8pc. 

Initially, the system is in thermal and near turbulent 
pressure balance, i.e for turbulent pressure balance 

{v\ = P -{v 2 ) w (1) 

pa 

holds, with velocities v and densities p. The indices stand for 
the cold and warm phase. The near pressure balance reduces 
motions due to pressure differences and therefore allows us 
to make more valid statements about turbulent transport or 
mixing in this idealized setup. 

The numerical scheme is based on the 2nd order 
Bhat na gar-Gross-Krook f or malis m ^jjPrend*rgast_&2£yl 
[ 19931: ISlvz fc PrendergasH Il999t iHeitsch et alJ 120041 : 
ISlTzetani2005|) . allowing control of viscosity and heat 
conduction. The code evolves the Navier-Stokes equations 
in their conservative form to second order in time and 
space. The hydrodynamical quantities are updated in time 
unsplit form. 

We employed the same heatin g and coo ling prescrip- 
tions as IHeitsch et all J20051 l200rj) , based on IWolfire et ail 
(1995). The same caveats apply, especially that we are dis- 
cussing the mixing between the warm and cold ISM, and 
that we are not including molecular gas. Thus, while the 
parameters of our cold cloud are consistent with values for 
"Giant Molecular Clouds", we neglect the effects of molec- 
ular line cooling and chemistry. 
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Figure 1. Column density maps of model A2b (top), C2b (cen- 
ter), and A21 (bottom) projected along the z-axis, for times 
t = T e , 3t<; and 5r e , where r e ~ 2Myrs is the initial turbulent 
crossing time. Color corresponds to log 7V[cm -2 ]. 

The code is equipped with Lagrangian tracer particles 
that are initially deployed within the cold cloud at a resolu- 
tion of one particle per grid cell. The particles are advected 
with the gas flow, so that they allow us to study the history 
of the cold gas. 

We restricted the models to hydrodynamics with heat- 
ing and cooling, leaving out gravity and magnetic fields. 
Depending on their strength, fields could suppress shear in- 
stabilities, while gravity might lead to more compact dense 
structures and fragmentation. This could have a twofold ef- 
fect, as will be discussed in ^1 



3 RESULTS 

3.1 Morphologies 

A first impression of the efficiency of turbulent mixing can 
be gleaned from a time sequence of column density maps 
(Fig. 0. Column densities are integrated along the z-axis 
and shown at times t = T e , t = 3r e and t — 5r e , where 
r e = L/vo ~ 2Myrs, a nominal turbulent flow crossing time 
in terms of the total box length and the rms velocity in 
the warm diffuse gas. Despite the fact that we are looking 
at column densities here, the overall effect of the cloud's 
turbulent dispersal is clearly visible. 

The most noticeable difference between the adiabatic 
models A2b and A21 (top and bottom) on the one hand 
and the radiative models (C2b, center) on the other is that 
for the latter, the transition between low and high column 
densities is much more marked, i.e. the column density maps 
appear (especially in the later stages) less "frothy". As we 
will see below, this is a direct consequence of the cooling. 

In overall appearance, for both the adiabatic and radia- 
tive models the cloud stays more or less compact, although 



fingers of cold, dense material are clearly sprouting from the 
core. 



3.2 Mass fractions and Dynamics of the Warm 
and Cold Gas 

Since we will be discussing the turbulent mixing of the 
warm and cold gas, first, we need to understand the dy- 
namics and evolution of each of the phases. Note that - 
strictly speaking - the concept of "phases" can be mislead- 
ing not only for the adiabatic case, but also for the thermally 
bist able case, because of the imp ortance of dynamics (see 
e.g.lVaz auez-Semadeni et al . 2000:: IVaza uez-Semadeni et alJ 
l200fJ: iHeitsch et alJl200fj) . Motivated by the bistable mod- 
els (sequence C), we split the model-ISM into three regimes, 
namely a cold phase with temperatures 300K > T, a warm 
phase (3000K < T) and an intermediate phase with 300K < 
T < 3000K. The corresponding mass fractions (Fig. |21 top 
panel) evolve quite differently for adiabatic and radiative 
models. The mass fractions are taken over the whole simu- 
lation volume. Since the cloud in model A21 has an initial 
radius twice as large as in the other models, it starts off with 
a larger cold mass fraction and a smaller warm mass frac- 
tion. The intermediate temperature regime evolves similarly 
to that in models A2a/b. 

In the adiabatic case (A2a/b/l), some of the cold gas is 
lost to the intermediate regime, while the mass fraction in 
the warm phase stays pretty much constant. In contrast, the 
radiative case (C2a/b) keeps constant mass fractions in each 
of the three temperature regimes over the whole simulation 
time (i.e. 5 dynamical times). Heating and cooling timescales 
are much shorter than the dynamical timescales, so that gas 
cannot collect in the intermediate regime, which for model 
C2a/b corresponds to the thermally unstable regime. In par- 
ticular, in C2a/b gas which is heated by (viscous) shear at 
the cloud boundaries "falls back" to its previous (cold) ther- 
mal state, while for A2a/b, this gas has no way to return to 
its previous temperature except by adiabatic cooling. Vary- 
ing the resolution does not change the mass fractions ap- 
preciably. Model C2a/b displays a tiny increase of the cold 
mass fraction. This stems from the compression of warm gas 
when it collides with the cold cloud rim. Once the warm gas' 
density increases, cooling sets in, and this compressed gas 
is added to the cold gas component. However, as Figure [5] 
shows, this effect is negligible. 

In contrast to the evolution of the mass fractions differ- 
ing for the adiabatic and radiative models, the rms velocities 
of the gas in the cold and warm temperature regimes evolve 
similarly for both sets of models. Figure [5] (bottom panel) 
mirrors the initial pressure balance (eq.0: velocities in the 
cold phase start out lower by a factor of approximately 17. 
For each of the models A and C, the velocities in the cold 
and warm phase decay, albeit at different rates. Resolution 
effects do not affect the decay (the thin and thick lines are 
nearly indistinguishable). Comparing model A to C, the ra- 
diative losses occurring in C do not lead to significantly dif- 
ferent decay rates. Radiative losses would become important 
in regions of high compression, however, the turbulence ini- 
tially decays quickly below Mach 1 in both the warm and 
cold gas, limiting the compression. We discuss the time evo- 
lution of the cold and warm pressure profiles in detail in 
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Figure 2. Top: Mass fraction over the whole simulation volume 
for the three thermal regimes indicated, against time, for models 
A2a/b and C2a/b at resolutions N = 128 3 and N = 256 3 , and for 
model A21 (at N = 256 3 ). Bottom: Rms velocity (v 2 ) 1 / 2 for cold 
and warm gas against time, for the same models and resolutions. 
Note that red and blue lines here denote T > 300K and T < 300K 
respectively. 



The larger cloud (model A21) shows the same velocity 
decay as its small counterparts. 

3.3 Radial Profiles 

Does the cloud break up in the turbulent environment, and 
if so, how quickly does this happen? The answer depends 
strongly on the quantity we are looking at. Since we are 
interested in an average measure of the cloud's structure, 
we take averages over shells and discuss the resulting radial 
profiles. The shell centers coincide with the instantaneous 
center of mass of the cloud (the cloud itself moves a little 
in the background flow). We begin with the density profiles 
(Fig.®. 

Clearly, the turbulent motions lead to a spreading-out 
of the cloud, mirrored in the density-weighted average radius 



(R)i 



/rf rp(r)dr 
/rf P(r)dr 



(2) 



of the cloud, which roughly doubles within 5r e = lOMyrs. As 
Figure already indicates, the radial density profiles show a 
substantial variation at fixed radius R (Fig El lower panel), 
up to six times the actual density value. The quantity 

a p (R)/p(R) = ((p(R) - {p)) 2 ) 1/2 I p(R) (3) 

plotted is the relative standard deviation on the mean den- 
sity over a sphere at radius R. The strongest variations are 
expected at the cloud rim, which travels (upper panel) out- 
wards, so that the peak of the density variations is seen 
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Figure 3. Top: Radial density profiles for models A2b, C2b and 
A21, at the times indicated. Model A2b is shown in thick lines 
in both plots. The thin solid step function gives the initial den- 
sity profile for models A2b and C2b. The initial density profile of 
model A21 extends out to 8.8pc. Bottom: Relative standard de- 
viation of density profile (eq. |3. The dotted vertical line denotes 
the initial radius. 



at larger radii for later times. Already after one dynamical 
time, the cloud is far from being a solid sphere. 

A more stringent measure is the radial volume filling 
factor profile for gas in the cold phase (Fig.^J. Volume filling 
factors are measured on shells at given R. As in >I3.2I the 
temperature threshold to distinguish between the warm and 
cold phase is set at T = 300K. The step function in the upper 
panel gives the initial condition, which of course shows a cold 
gas filling factor of 1 for R ^ 4.4 pc (or R $J 8.8 pc for model 
A21). 

The first obvious difference between the radiative and 
adiabatic case is that the adiabatic case seems to mix faster 
the cold and warm phase at more radii. However, this is not 
that surprising, since for the adiabatic case any gas with 
T > 300K is lost for the cold gas filling factor, while for the 
radiative case, gas which has left the cold phase can only be 
found in the warm phase (there is no intermediate-phase gas, 
Fig- HI top panel) but because of the short cooling times, this 
gas can quickly return to the cold phase. In other words the 
filling factor gives an unambiguous measure of the degree of 
mixing between cold and warm gas only for the thermally 
bistable, radiative case (C2b), since there is (close to) zero 
conversion between the gas phases (see also discussion on 
tracer particles below). 

After 2 Myrs (corresponding to t = r e ), the volume fill- 
ing factor of the cold gas measures 40% for the radiative 
case, and ~ 25% for the adiabatic case at the initial cloud 
boundary, i.e. more than half of the volume is occupied by 
warm gas. Note that especially in the radiative case, the 
mass fractions of cold and warm gas stay constant, i.e. cold 
and warm gas are bodily transported. This can be gleaned 
from Figure |^| It shows the histogram of the fraction of 
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Figure 5. Histogram of fraction of tracer particles in the flow 
staying at T < 300K for a time interval At. Most tracer particles 
stay cold over the whole simulation time (At > 9Myr) for model 
C2b (thick line), while for model A2b (thin line), only ~ 60% of 
the particles stay at T < 300K. 



Figure 4. Top: Radial cold gas filling factor profile for models 
A2b, C2b and A21 at the times indicated. Model A2b is shown 
in thick lines in both plots. The thin solid step function gives 
the initial cold gas filling factor profile for models A2b and C2b. 
The initial profile of model A21 extends out to 8.8pc. Bottom: 
Relative standard deviation of filling factor profile (analogous to 
cq.[3J. The dotted vertical line denotes the initial radius. 

Lagrangian tracer particles Np(T < 300K)/7Vp within the 
cloud that stay at temperatures T < 300K for a time inter- 
val At. If there were no turbulence in the models, all the 
particles would stay at T < 300K for the whole duration of 
the simulation, i.e. we would have Np(T < 300K)/Np = 1 
at At = 10. 

Most of the particles stay cold for longer than 9Myr, 
i.e. for the simulation's time extent. In other words, the con- 
stancy of the mass fractions in Figure |5] indeed results only 
to a very minor extent from the conversion of cold to warm 
gas and vice versa: mass in the cold and warm regime is sep- 
arately conserved. This is less valid for the adiabatic case, 
which obviously loses some of its cold material. But still, 
the bulk of the initially cold gas stays cold for the whole 
simulation time. 

Besides turbulent transport, an overpressure in the cold 
cloud could lead to "expansion" . So far, our diagnostics can- 
not distinguish between these two mechanisms for cloud ex- 
pansion. Figure |S| shows the radial profiles of the pressures 
for models A2b, C2b and A21. The total pressure Ptot has 
been split into the thermal pressure P, and the turbulent 
pressure Ptrb, from which we have removed the contribution 
of the translational velocity (v), as 

Ptot = P + p(v-(v}) 2 , (4) 

where the average extends over coherent cold and warm re- 
gions. The top row gives the radial pressure profile for gas 
with T > 300K, and the bottom row shows the profile for 
gas with T < 300K. The time sequence reveals the mixing of 
the warm and cold component, since P(T > 300K) can be 
defined at smaller radii for later times, while P(T < 300K) 
spreads outward to large radii resulting in warm and cold 



gas co-existing at an increasingly larger radial range. For 
t = 0, models A2b and C2b are of course identical. The 
thermal pressure is constant, while the turbulent pressure 
(and thus the total pressure) is slightly lower within the 
cold cloud than in the warm medium. This serves as a safe- 
guard against initially overpressuring the cloud by turbulent 
pressure. At t — r e (center column), the cold dense material 
(i.e. the cloud) in both models is overpressured relative to 
the warm gas. This comes mostly from an overshoot in the 
density, since turbulent pressure and thermal pressure have 
(approximately) the same radial dependence: the cloud is 
initially slightly compressed by the higher turbulence in the 
warm medium (see also the radial density profiles, Fig. |3J. 
In Figure (bottom panel) we saw that the rms velocity - 
and because of the close to constant density in the respec- 
tive phases, the kinetic energy - decays faster for the warm 
gas than for the cold gas. Both effects together lead to the 
overpressure in the cold gas at t = T e . At later times, this 
pressure imbalance has canceled out, and the warm and cold 
phases are mixed. There is still a turbulent pressure excess 
at small radii, however, the turbulence has decayed so far 
that its dynamical effects are insignificant. 

For the adiabatic run, at t = T e the thermal pressure 
has risen from its initial value both for the cold and warm 
gas. The turbulent pressure on the other hand, has dropped 
by about a factor 3 for the warm gas, but not as much for 
the cold gas, suggesting that turbulent energy in the warm 
phase has been used to heat the warm phase and to drive 
the turbulence in the cold phase, since the thermal energy in 
the cold phase increased and the turbulent energy remains 
unchanged. With the total energy conserved, some of the 
energy in the cold gas has to come from the warm gas. For 
the radiative model at t = r e the transfer of energy from 
the warm to the cold gas is less marked but seems to have 
occurred nevertheless. The cold gas' thermal and turbulent 
pressure in the innermost radii increase from their initial 
values. The thermal pressure of the warm gas on the other 
hand has hardly changed from its initial value even though 
the turbulent pressure has dropped by almost an order of 
magnitude from its initial value at all radii where warm 
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Figure 6. Radial pressure profiles for models A2b (top), C2b 
(center) and A21 (bottom) for times 0, T e and 5r e . In each panel, 
from top to bottom: radial pressure profile for gas at T > 300K, 
total radial pressure profile, and radial pressure profile for gas 
at T < 300K. Plotted are the total pressure (solid lines), the 
thermal pressure (dashed) and the turbulent pressure (see text, 
dash-dot). The three columns per model denote the measurement 
times t = 0, t = T e and t = 5t.. 
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Figure 7. Left: Radial optical depth profile for models C2b (thin 
lines) and A2b (thick lines), at the times as indicated in the plot. 
The thick dotted line gives the initial optical depth profile, while 
the thin dotted line denotes t v = 1. Right: Same plot for models 
A21 (thin lines) and A2b (thick lines). Note that the optical depth 
range is larger for model A21. 



gas exists. Again, the thermal and/or turbulent pressure in- 
crease in the cold gas has to come from the warm gas. Model 
A21 shows essentially the same behavior as model A2b. 

The radial density profile allows us to determine the av- 
erage optical depth at a given radius (Fig. 0. To arrive at 
meaningful optical depths, we scaled the effective absorption 
coefficient ct v such that the central optical depth is initially 
arbitrarily set to r c = 7.5 for models A2a/b and C2a/b, and 
to r c = 15 for model A21. This corresponds to an absorption 
coefficient of a v — 3.6 x 10~ 21 cm 2 . For comparison, Spitzer 
(1961) gives an effective combined cross section for absorp- 
tion and scattering on grains of a v = 1.2 x 1CP 21 cm 2 at a 
wavelength of 100 nm. Likewise, we do not take into account 
the effect of scattered light. 

There are no marked differences between the adiabatic 
and the radiative case. This as well as the fact that R(r = 1) 
does not change appreciably over time is a direct conse- 
quence of the radial density profiles (Fig. [3J. At t = r e , the 
optical depth in the center has increased due to an initial 
compression, just to drop at later times because of turbu- 
lent dispersal. This "overshoot" is mirrored in the density 
profiles (Fig. [3] top panel) . 

The central optical depth r c drops by approximately a 
factor of 2, but does not fall below r = 1. That r c changes at 
all with time might come as a surprise, but is a consequence 
of the exchange of dense and diffuse material at approximate 
pressure equilibrium: the filling factor of dense material on 
a spherical surface at fixed radius decreases, lowering the 
central optical depth. 

3.4 Turbulent Diffusivity 

The evolving radial density and filling factor profiles (Figs.|3] 
and^J suggest to estimate the (turbulent) diffusivity by fit- 
ting them to the expected profiles resulting from the time 
evolution of a step function under the effect of diffusion. 
An (inert) quantity q in a turbulent environment obeys the 
advection-diffusion equation 

(d t + u ■ V)g = AV 2 g, (5) 

where u is the velocity, and A the microscopic diffusivity, 
which has been assumed to be independent of location and 
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direction. Equation Elcan be rewritten (see e.g. appendix in 
iHeitsch et afl l2004) as 



d t (g) = (A e +A)V 2 < 9 ), 



(6) 



under the - contestable - assumption of a separation be- 
tween the small scale turbulent velocity field and the large 
scale variations in the quantity (q) , where the averages have 
removed variations due to the small-scale turbulence. The 
"turbulent diffusivity" A e = u rms L is the product of the 
rms velocity and the characteristi c length scale over which a 
gas parcel maintains u rma (see e.g. lLandau fc Lifshit J l966). 
In the ISM, the Reynolds number Re = A e /A ^> 1, generally, 
so that A can be neglected for turbule nt transport s tudies. 
Quasi-linear diffusion theory (see e.g. iMoifal ll 8 ) holds 
that A e can be regarded as constant. 

Since our model clouds start o ut with a uniform den- 
sity, w e can follow the discussion bv lde Avillez fc Mac Lowl 
(2002) and study the turbulent diffusive evolution of a 
step function profile. A one-dimensional density distribution 
evolves as 

n {x, t) = t. / n{x',t = 0) e~ (x - x > /AXct dx (7) 

under diffusion. The initial conditions can be written as 
nix Ro,t = 0) = ni and n(x > Ro,t = 0) = no, where 
Ro is the initial cloud radius at t = 0. Then, the density 
distribution at time t is given by 



i(x,t) = n + 



(8) 



where erf(i) is the error function. A similar expression is 
valid for the filling factor f(x,t), if we replace n(x,t) by 
f(x, t), and set no = and m = 1. Since the initial radius _Ro 
is known, we can fit equation[S]to the available density and 
filling factor profiles at given times t, and thus determine 
the diffusivity A e by a I-para meter Levenberg-M arquardt 
least-squares minimization (e.g. IPress et al.lll992T) for each 
available time step of a model. 

Figure |S] shows the resulting diffusivities for the filling 
factor profiles (top left) and the density profiles (top right) 
for all models. The center row gives the l-a errors of the 
profile fits, and the bottom row contains the reduced Xred 
for the filling factor profiles (bottom left) and for the density 
profiles (bottom right). Reliable fits we select (somewhat 
arbitrarily) by Xred < 2, agreeing well with a selection by 
eye. 

Clearly, the filling factor profiles lead to much better 
fits than the density profiles. Figures |3] and 2] explain this: 
the filling factor cannot rise above 1 or drop below 0, thus 
constraining the profiles for the least-square minimization, 
while the density profiles actually increase above ni and 
drop below no because of traveling waves. Moreover, in de- 
riving equation [H] we assumed that n(x ^ Ro,t = 0) = m 
is valid for all x < Ro, and not only for ^ x ^ Ro, as 
by construction of the initial conditions. This assumption is 
certainly not valid any more at later times. 

We first notice that the diffusivities determined by least- 
square fitting (symbols in top row of Fig- |HJ are n °t constant 
with time. Until ~ 3Myrs they all increase. The increase 
however is less than an order of magnitude. The filling factor 
profiles return slightly decreasing dif fusivities for later times 
and for all models except model A21. lde Avillez fc Mac Lowl 



(2002) found that the diffusivities increase exponentially 
with time. Since they start with an unperturbed medium 
and then drive the turbulence via supernova explosions, the 
growing diffusivities could be a result of the increasing rms 
velocity. In our models, a similar effect is causing the ini- 
tial rise of A e . The turbulent diffusivity increases while the 
turbulent cascade is building up. At later times, the overall 
decay of the turbulent rms velocity then leads to a decrease 
of the turbulent diffusivity. 

The density profiles do not provide as nearly a complete 
picture as the filling factor profiles. Reliable fits are difficult 
to get because of shock waves traveling through the cloud. 
The few reliable fits clearly show that the diffusivities are 
not constant, but range within the same values as the ones 
derived from the filling factor profiles. As a crude check, the 
classically expected diffusivities 



(R) 



M U r 



(9) 



are shown in the top right panel of Figure [8] in lines. Ini- 
tially, they decrease slightly, mirroring the decay of u rma , 
while the growing cloud radius compensates for the decay- 
ing velocity at later times. The values are within a factor 
of 2 at most times with the diffusivities derived from the 
profile fits. A X e — 3 x I0 23 cm 2 s -1 would correspond to 
turbulent transport at I km s _1 over I pc. 

Clearly, the dif fusivities are not constant, which a grees 
with the findings oflde Avillez fc Mac Lowl <l2002l . l2003ft and 
iKlessen fc Linl ll2003l) . In a fully turbulent medium - whose 
energy scale distribution follows at least qualitatively a tur- 
bulent spectrum - the scale separation between turbulence 
and quantity to be diffused does not exist. 



3.5 Cloud Brightness 

Turbulence considerably changes the density profiles of our 
model cloud, and reduces the optical depth at the center 
of the cloud. Does turbulent transport ( H3.4H and the ex- 
change of warm diffuse and cold dense material ( H3.3t carve 
tunnels and holes in the cloud through which radiation can 
enter? For an answer, we have to de t ermin e the radiation 
field within the cloud. iBethell et al.l (|2004) discussed this 
question with the help of a spherical cloud inscribed in 
a periodic box of evo lved self-gravitating MHD-turbulence 
dHeitsch et alj|200fal lbT). This allowed them to study the ra- 
diation field inside an evolved structured cloud. Here, we 
are interested in the timescales on which the radiation field 
changes. 

To get a measure of the brightness inside the cloud, we 
determine the intens ity of the radiation field at each point 
inside the cloud fsee lBethell et aDl2004l) . Figure |5] shows a 
scatter plot of the intensities inside the volume occupied by 
the original cloud against radius, for the four times indi- 
cated. 

The brightness is calculated for each point by measuring 
the incident radiation for a given number of rays and aver- 
aging over the resulting sky. The ray number is determined 
such that at the outer surface of the cloud each resolution 
element is hit by one ray. Note that while we plot only points 
up to the original radius of the cloud, the radiative trans- 
fer includes the whole domain, i.e. we do not lose "dense" 
material outside R = 4.4pc which could shadow the inner 
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Figure 8. Diffusivities (top row), l-cr errors for profile fits (center row), and reduced \ 2 f° r profile fits (bottom row), against time for all 
models. Left panels show results derived from filling factor profiles and right panels show results derived from density profiles. Symbols 
in the top row panels denote diffusivities derived from profile fits (eq. |SJ, lines in the top right hand panel stand for A e derived from 
cq. [5] Only fit diffusivities for Xred ^ are snown ' 



parts. For extinction, we consider an effective extinction co- 
efficient. Actual scattering is not included. In that sense our 
brightness estimates are lower limits. 

The radially binned intensities resulting from plots like 
Figure [5] for models A2b, C2b and A21 are shown in Fig- 
ure ED Thick lines denote spherical volume averages at con- 
stant R, and thin lines spherical mass averages. Without 
self-gravity, the density contrasts in the cold material are 
small, so that volume and mass averages do not differ widely. 
Clearly, Figure I1UI mirrors the effect observed in Figure |5] 
the cloud gets "bright" within a few dynamical times, i.e. 
the turbulence opens holes. Even for the larger cloud (A21), 
the relative intensity does not drop below ICF 3 anywhere in 
the cloud. 

At late times, even the innermost regions receive more 
than a tenth of the incident radiation (models A2a/b, 
C2a/b, Fig. lllH : turbulence (indirectly) lights up the cloud. 
The intensity at the center of the cloud increases by at 
least 2 orders of magnitude. The central intensity for the 
models at lower resolution (thin lines, A2a and C2a) grows 
more smoothly than the intensity at higher resolution, and it 



reaches a slightly higher value at the end of the simulation. 
These are both resolution effects: the larger scatter comes 
from a more structured velocity field, and the larger central 
intensity is a consequence of fewer available grid cells along 
the line of sight. Still, for model A21 - with its initial cloud 
radius twice the size of those in models A2a/b -, the central 
intensity increases by four orders of magnitude (compared to 
two for models A2a/b and C2a/b). Thus, we do not expect 
higher resolution to lead to progressively smaller central in- 
tensities. 

The central intensities discussed so far are derived from 
averaging the incoming radiation over the whole "sky" as 
seen from the cloud center. If the turbulence digs tunnels in 
the cloud through which the irradiation can enter deep into 
the cloud, then this should be mirrored in the minimum (and 
maximum) optical depth at the center of the cloud (Fig. ll2ll . 
The extrema of the optical depths were taken over the whole 
sky as seen at the center of the cloud. The mean optical 
depth corresponds to the center intensities discussed above 
(note however that in order to derive the center intensities, 
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Figure 10. Radially binned intensities (see EI3.5I for models A2b (left), C2b (center) and A21 (right). Line styles denote the evolution 
time as indicated in the diagrams. Error bars represent errors on the mean. Thick lines arc volume-averaged intensities, thin lines are 
mass-averaged intensities. 
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Figure 12. Optical depth at the center of the cloud against time. Left column: Minimum optical depth. Center: Maximum optical depth. 
Right: Average optical depth. Top row: Models at resolution 128 3 (thin lines) and 256 3 (thick lines). Bottom row: Parameter study. Thin 
lines denote models with an initial rms velocity of Mach 2, while thick lines stand for those with Mach 3. 



the intensities are averaged over the sky, not the optical 
depths). 

As with the intensities (Fig. II It , the lower resolution 
runs lead to smaller average optical depths (Fig. 1121 top 
right). This is mostly a consequence of the lower minimum 
optical depths (top left of same figure) and of the dominance 
of large-scale motions. Thus, the models run at 128 3 grid 
cells are not fully resolved, but should be used for demon- 
stration purposes only. The strong spike for model C2b is 
a direct effect of the radiative cooling: compressed regions 
during the early turbulent evolution can reach much higher 
densities than possible for the adiabatic case. Although the 
minimum and maximum optical depth starts out at the 
same value, the maximum optical depth increases due to 
local compressions. Of course, the directions of minimum 
and maximum optical depth change with time. Models with 
different initial Mach number vary less than those with same 
Mach number but run at different resolution. Since the tur- 



bulence in our models is decaying, this is not surprising: 
most of the energy is lost in the early stages of the evolu- 
tion. 

Summarizing Figure 1121 the minimum optical depths 
drop by a factor of 4 to 8 (depending on the model), while 
the maximum optical depths increase by a factor of 2. Thus, 
while most of the brightening of the cloud stems from tur- 
bulent transport of dense material (see H3.2l and 13.311 . com- 
pression and corresponding evacuation contributes as well. 
Finally, the average optical depths drop by a factor of 4 or 
more. 



4 SUMMARY 

Motivated by the highly dynamical nature of molecular 
clouds and the question of how long a cold cloud can survive 
as a well-defined entity in a turbulent environment, we in- 
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Figure 9. Scatterplots of relative brightness within cloud A2b, 
against radius, and for times as indicated in panels. 
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Figure 11. Central relative intensity against time for models 
A2a/b, C2a/b and A21. Because of the larger cloud radius, model 
A21 stays at lower intensities. 



vestigated the evolution of a cold cloud in a warm, turbulent 
medium, assessing the efficiency of the mixing by measuring 
the brightness inside the cloud. We showed that the proper- 
ties of the initially solid cloud change thoroughly within a 
few dynamical timescales. 

While the column density maps can mislead, the bright- 
ness distribution inside the cloud clearly demonstrates that 
within a few dynamical times, the cloud becomes completely 
porous. Specifically, the cloud radius doubles within 5 dy- 
namical times, while the filling factor of the cold gas drops to 
less than 50% at the original cloud radius within one dynam- 
ical time. Cold material is bodily transported and exchanged 
with warm material (Figs. [5] and Q)J. The turbulent diffusiv- 



ity A e (eq. [S], [5] and Fig. |HJ is consistent with a fiducial 
number of X e « 10 23 cm 2 s— 1, corresponding to turbulent 
velocities and length scales of approximately 1km s _1 and 1 
pc. The diffusivities are not constant with time. 

The brightness estimates are lower limits, since we do 
not include scattering in the radiative transfer. Gravity 
might change the results by leading to additional fragmen- 
tation and thus to a growing discrepa ncy between volume - 
weighted and mass- weighted intensity tBethell et aljl2004f) . 
However, volumewise, the cloud would get even brighter this 
way due to the smaller filling factor of the dense gas. The 
effects of H2 (re-)formation remain to be discussed in sub- 
sequent models. If the turbulence were continually driven, 
we would expect an even faster dispersal. In that sense also, 
the presented timescales are only upper limits. 

The high radiation field within the cloud could strongly 
affect the chemistry and the dynamical state of the cloud 
ide Boisanger fc Chiezelll99ll) . leading to additional heat- 
ing, ^-destruction, and thus faster cloud disruption. 
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